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Abstract 

We present a short summary of parton saturation concepts as seen in 
deep inelastic scattering. 

1 Introduction 

The deep inelastic scattering (DIS) experiments, in which leptons probe nucleons with the help 
of electroweak bosons, reveal that nucleons consist of partons. These are colored quarks of 
Quantum Chromodynamics (QCD) which carry approximately half of the nucleon's momentum. 
The missing half is provided by gluons to which the electroweak bosons do not couple. Thus, 
although not directly probed, gluons are extremally important for the description of the nucleon 
structure. Quantitatively, this is summarized by the DGLAP evolution equations of QCD which 
govern the dependence of the quark and gluon distributions in a nucleon on a scale (identified 
in DIS with photon's virtuality = —Q^). The sign of the logarithmic derivative, dF2 / d log Q^, 
at different values of the Bjorken variable x is determined by the relative contribution of quarks 
to gluons. In the limit x — 0, studied intensively by the experiments at HERA, the deep inelastic 
processes are dominated by a strongly rising gluon distribution. Therefore, in the small-x limit, 
gluonic systems inside the nucleon are predominantly studied. The description of processes in 
such systems, using perturbative QCD (pQCD), is the aim of this presentation. 

2 CoUinear factorization versus /ct -factorization 

In the electron-proton DIS, the measured proton structure functions. Ft and Fl, are related to 
the parton distributions through the colUnear factorization formula resulting from pQCD: 

A^"^ ix a ) 

i=q,q,g n=l ^ 

where (g) indicates integral convolution in parton longitudinal momentum fractions, as = as{Q'^) 
is the running strong coupling constant, c!^^^{z, Os) are perturbatively computed coefficient 
functions and fi{x,Q'^) are quark, antiquark and gluon distributions (multiplied by x). The 
Q^-dependence of the parton distributions is determind by the DGLAP evolution equations [1] 
with initial conditions which are fitted to data. The first term on the r.h.s. of eq. ^ provides the 
leading twist-2 description with logarithmic dependence on while the remaining terms, called 
higher twists, seems to be suppresses for large Q^. In the standard analysis, a global fit of the 
leading twist formula to the HERA data on F2 = Ft + Fl, together with cross sections of other 
hard processes, leads to the determination of the parton distributions shown in Fig. [T] A distinct 
feature of this determination is a strong rise of the gluon and sea quark distributions for x ^ 0. 
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Fig. 1: Parton distributions from a global fit to the HERA data as functions of x for fixed — 10 GeV^. 

A closer theoretical examination of the small-x scattering reveals that for not too high Q^, 
the higher twist terms cannot be neglected since they are enhanced by powers of as log(l/2;), 
when the smallness of Ug is compensated by large logarithm of x. The relevant resummation of 
such terms in the leading (LO) and next-to-leading (NLO) logarithmic approximation leads to 
the BFKL approach to the structure functions with the following /c^-factorized form [2]: 



where the impact factor ^{kj^/Q'^,as{kT)) describes the interaction of the virtual photon with a 
gluon with nonzero transverse momentum kx- In the LO this is the process: 7* {Q'^)g{kT) — > qq- 
The function f{x,kT) is called unintegrated gluon distribution which obeys the BFKL equation 
[3] and is related to the gluon distribution g{x,Q'^) through the formula 



From the solution of the BFKL equation, the small-x limit is dominated by the gluon distribution 
with the power-like rise, /(x,A;t) ~ and A ^ 0.3. There is a general agreement, based 
on the experience with the Froissart-Martin bound, that such a rise of the gluon distribution, and 
in consequence F2, violates unitarity and eventually must be tamed. The BFKL solution is also 
plagued by diffusion to infrared, namely, the -integration in the pQCD formula Q is quickly 
dominated by the contribution from the soft momenta region, kr ~ ^qcd, where the Landau 
pole of asikx) is encountered. A cure for these problems is absolutely necessary. 

3 Parton saturation 

The taming of the power-like rise of the gluon distribution xg{x, Q^) was addressed for the first 
time by Gribov, Levin and Ryskin in [4] in the double logarithmic approximation. Summing fan 
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Fig. 2: Saturation line in the {x, Q^)-plane. 



diagrams, which take into account the fusion of i-channel gluons, the Unear DGLAP equation 
for the gluon distribution receives a negative, nonlinear term. 



ain(l/rE)ainQ2 



vr- 



(4) 



where = NcOs/tt and the parameter R controls the strength of the nonlinearity. With such a 
modification, the gluon distribution saturates for x — > 0, and so does the structure function. This 
result was extended in [5] by including nonlinear modifications for the sea quark distributions. A 
crucial feature introduced by the nonlinearity is an x-dependent saturation scale Q^{x), defined 
as a value of Q"^ for which the nonlinear term in eq. ^ is comparable with the linear one: 
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Therefore, saturation effects are important when the number of gluons per unit of rapidity, xg, 
times the gluon-gluon interaction cross section, as/Q"^, approaches the geometric size of the 
nucleon or a gluonic system inside the nucleon ("hot spot")- In such a case, a simple additive 
treatment of parton emission breaks down and gluons start to annihilate. Since from ^ xg 
and xg ~ x^'^ before the saturation limit is reached, we find that » ^qcd sufficiently 
small X, and the presented approach based on perturbative QCD is justified . This is schematically 
illustrated in Fig. |2] where two regions separated by the saturation line, = Ql{x), are shown. 
Below this line, in the dilute region, the linear evolution equations are valid, while approaching 
the line, the saturation region is entered with nonlinear equations describing parton saturation. 

Eq. (|4l) is a rather crude approximation since it is valid in the extreme case, x — > and 
oo. In the fcy -factorization approach the latter limit is relaxed and only large logarithms 
log(l/x) are relevant. Summing BFKL pomeron fan diagrams with triple pomeron vertices in 
the leading logarithmic approximation and in the limit of large number of colors Nc, the Balitsky- 



Kovchegov (BK) equation for the unintegrated gluon density (/)(x, kx) is foun({3 [6,7]: 

dY(t){x,kT) = asx{-dL)4' - 0^ (6) 

where Y = log(l/3;) is rapidity, L = log A:^ and x is the BFKL characteristic function. This 
nonlinear equation generalizes the linear BFKL equation. The properties of its solutions were 
intensively studied both analytically [8] and numerically [9, 10]. The most fruitful approach is 
based on the relation to the known from statistical physics Fisher-Kolmogorov equation, which 
admits travelling wave solution. In our language, it means that the BK solution develops a satu- 
ration scale, Qs{x) ~ with known value of A [8], such that for small x we have 

(l)ix,kT) = HkT/Qsix)). (7) 

This property, called geometric scaling, was observed in the data from HERA [11]. Looking 
more carefully, for kx ^ Qs{x) the gluon distribution <j) ~ l/kj., while for small transverse 
momenta, Aqcd kr < Qs{x), the behaviour changes to logarithmic, ~ ln{Qs{x)/kT)- 
This is the illustration of the transition to saturation, when both the power-like growth in x and 
infrared diffusion in kx of the gluon distribution are tamed [10], see Fig. |2] with = kj^. 



4 Color dipole approach and beyond 

A more intuitive approach to parton saturation is provided by the color dipole approach [12, 13]. 
In the target rest frame, the DIS at small x can be formulated as the eikonal scattering of a color 
quark-antiquark dipole, formed by the splitting 7* — s- qq, on the target color field. The dipole 
scattering amplitude A^(x, y) is given by two Wilson lines coUinear to quarks velocity u 

iV(x,y) = l-^Tr[/(x)C/t(y), U{x) = F exp f^ig ^ dXu ■ AiXu + x)^ (8) 

where x and y are two dimensional vectors of quark transverse positions, conserved during the 
collision, and ^4 is a target color field. The deviation of the classcal quark trajectory from the 
light-like line defines the change of N with rapidity Y, which leads to the new BK equation for 
the dipole scattering amplitude [6]. Its solutions fulfil the unitarity bound, < 1. When the 
dependence on the impact parameter, b = (x + y)/2 is neglected, the new equation is equivalent 
to eq. (O after Fourier transforming of N/r"^ with respect to r = x — y. The BK equation in the 
transverse space was also obtained in the Mueller's dipole approach [13] in which the qq dipole 
develops a system of dipoles (by radiating soft gluons in the large Nc approximation) which 
subsequently multiple interact with a large nucleus target [7]. 

The dipole scattering amplitude is the basic ingredient in the computation of the nucleon 
structure functions at small x. In the last ten years, this amplitude was also modelled using the 
properties of the BK solutions such as color transparency, A^ ~ for a small dipole size r = |r|; 
geometric scaling, A^ = N{rQs{x)); and the unitarity bound, A^ < 1. A recent comprehensive 
review on the dipole models of DIS processes is presented in [14]. 

The BK equation describes unitarity corrections in the asymmetric configuration when 
the target is extended and dense and the projectile is small and dilute. In a more symmetric 

'(j> is related to the unintegrated gluon density / from Section|2]by f{x, k) ^ k^V\<f){x, k). 



configuration, e.g. in the pp scattering at LHC, the BK equation is no longer sufficient, which 
means that in the diagrammatic approach closed pomeron loops have to be taken into account 
besides fan diagrams. An interesting attempt in this direction was made in [15] where pomeron 
loops were modelled as color reconnections in the dipole cascades. The resulting scattering 
amplitudes respect the target-projectile symmetry and describes reasonable well the existing total 
and diffractive cross sections in the pp scattering. The pomeron loops were also studied in a 
statistical approach, based on the stochastic Fisher-Kolmogorov equation, finding a new kind of 
scaling called diffusive scaling [16]. Recently, high energy factorization theorems for the gluon 
production in heavy nucleus collisions were proven in the color glass condensate approach [17]. 

Acknowledgements This work has been supported by the Polish grant no. N N202 249235 and 
by the Research Training Network HEPTools (MRTN-2006-CT-035505). 

References 

[1] V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972); 

G. Altarelli and G. Parisi, Nucl. Phys. B126, 298 (1977); 

Y. L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977). 
[2] S. Catani, M. Ciafaloni, and F. Hautmann, Phys. Lett. B242, 97 (1990); 

S. Catani, M. Ciafaloni, and F. Hautmann, Nucl. Phys. B366, 135 (1991); 

J. C. Collins and R. K. EUis, Nucl. Phys. B360, 3 (1991). 
[3] V. S. Fadin, E. A. Kuraev, and L. N. Lipatov, Phys. Lett. B60, 50 (1975); 

L. N. Lipatov, Sov. J. Nucl. Phys. 23, 338 (1976); 

E. A. Kuraev, L. N. Lipatov, and V. S. Fadin, Sov. Phys. JETP 45, 199 (1977); 
1. 1. BaUtsky and L. N. Lipatov, Sov. J. Nucl. Phys. 28, 822 (1978). 
[4] L. V. Gribov, E. M. Levin, and M. G. Ryskin, Phys. Kept. 100, 1 (1983). 
[5] A. H. Mueller and J.-w. Qiu, Nucl. Phys. B268, 427 (1986). 
[6] I. Balitsky, Nucl. Phys. B463, 99 (1996). hep-p h/ 950 9348] 
[7] Y. V. Kovchegov, Phys. Rev. D60, 034008 (1999). | hep-ph/ 9901281[ 
Y. V. Kovchegov, Phys. Rev. D61, 074018 (2000). hep-ph/ 99 05214| 
[8] S. Munier and R. B. Peschanski, Phys. Rev. Lett. 91, 232001 (2003). |hep-ph/0309177[ 
S. Munier and R. B. Peschanski, Phys. Rev. D69, 034008 (2004). |hep-ph/0310357| 
S. Munier and R. B. Peschanski, Phys. Rev. D70, 077503 (2004). |hep-ph/0401215| 
[9] N. Armesto and M. A. Braun, Eur. Phys. J. €20, 517 (2001). hep-ph/ 104 03_8J 
[10] K. J. Golec-Biernat, L. Motyka, and A. M. Stasto, Phys. Rev. D65, 074037 (2002). |hep-ph/0110325| 
[11] A. M. Stasto, K. Golec-Biernat, and J. Kwiecinski, Phys. Rev Lett. 86, 596 (2001). |hep-ph/0007192| 
[12] N. N. Nikolaev and B. G. Zakharov, Z. Phys. €49, 607 (1991). 
[13] A. H. Mueller, Nucl. Phys. B415, 373 (1994); 

A. H. Mueller and B. Patel, Nucl. Phys. B425, 471 (1994). |hep-ph/ 940325 6| 
A. H. Mueller, Nucl. Phys. B437, 107 (1995). he p-ph/9408245| 
[14] L. Motyka, K. Golec-Biernat, and G. Watt (2008). 10 809 . 4191] 

[15] E. Avsar, G. Gustafson, and L. Lonnblad, JHEP 07, 062 (2005). |hep-ph/0503181| 
E. Avsar, G. Gustafson, and L. Lonnblad, JHEP 01, 012 (2007). •hep -ph/0610157| 
E. Avsar, G. Gustafson, and L. Lonnblad, JHEP 12, 012 (2007). 07 09. 13681 

[16] A. H. Mueller and A. I. Shoshi, Nucl. Phys. B692, 175 (2004). |hep-ph/0402193| 
S. Munier, Nucl. Phys. A755, 622 (2005). hep-ph7050 114 9 ; 

E. lancu, A. H. Mueller, and S. Munier, Phys. Lett. B606, 342 (2005). ihep-ph/0410018| 
[17] F. Gelis, T. Lappi, and R. Venugopalan (2008). 804.26301 

F. Gelis, T. Lappi, and R. Venugopalan (2008). 807 . 1306] 



